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Abstract 

We derive a factorization theorem for production of an arbitrary number of color- 
singlet particles accompanied by a fixed number of jets at the LHC. The jets are 
defined with the standard anti-Zc-r algorithm, and the fixed number of jets is obtained 
by imposing a veto on additional radiation in the final state. The formalism presented 
here is useful for current Higgs boson analyses using exclusive jet bins, and for other 
studies using a similar strategy. The derivation uses the soft-collinear effective theory 
and assumes that the transverse momenta of the hard jets are larger than the veto 
scale. We resum the large Sudakov logarithms a" log^""™" (pj./p^^*") up to the next-to- 
leading-logarithmic accuracy, and present numerical results for Higgs boson production 
in association with a jet at the LHC. We comment on the experimentally-interesting 
parameter region in which we expect our factorization formula to hold. 



1 Introduction 



Accurate predictions for processes with a fixed number of final-state jets are crucial for many 
LHC applications. A well-known example is that of a Higgs boson decaying to M^-bosons 
at the LHC PQ |2]. The background composition to this signal changes as a function of jet 
multiplicity. In the zero-jet bin the background is dominated by continuum WW production, 
while in the one-jet and two-jet bins, top-pair production becomes increasingly important. 
The optimization of this search requires cuts dependent on the number of jets observed, and 
therefore also on theoretical predictions for exclusive jet multiplicities. 

Theoretical predictions for processes with an exclusive number of jets are notoriously 
difficult to obtain. Fixed-order perturbation theory is plagued by large logarithms of the form 
\n{Q / pl^^") , where Q denotes the hard scale in the process, such as mn- For experimentally 
relevant values ~ 25 — 30 GeV, residual scale variations in fixed-order calculations lead 
to estimated errors that do not accurately refiect uncalculated higher-order corrections [U 
m [5] . Progress in resummation of these logarithms to all orders has been slow. Event-shape 
variables such as jettiness [0] allow resummation of jet-veto effects to arbitrary logarithmic 
accuracy, and have been applied to study the production of vector bosons or Higgs bosons 
plus multiple jets at the LHC [HE]. However, experimental measurements typically utilize 
jet algorithms such as the anti-fcr algorithm, and conclusions drawn from calculations using 
jettiness necessarily remain qualitative only. Resummation of jet-veto logarithms for the 
Higgs cross section in the zero-jet bin in the presence of the anti-fc^ algorithm has been 
performed at next-to-leading logarithmic (NLL) accuracy using the semi-numerical program 
CAESER [_5_£| Recent work has extended these results to their NNLL accuracy P HO] . It has 
been pointed out the potentially large In R corrections, where R is the jet-radius parameter 
in the anti-fc^ algorithm, could have a significant numerical impact on the predictions [11]. 
These terms have yet to be studied at all orders and warrant further investigation. 

We consider in this manuscript the resummation of the jet-veto logarithms for production 
of one or more color-neutral particles, such as a Higgs boson or electroweak gauge bosons, in 
association with one or more jets. We accomplish this by deriving a factorization theorem 
using soft-coUinear effective theory (SCET) [121 ttSl ttH ttSl [16] that assumes that the trans- 
verse momenta of the hard jets are larger than the veto scale. As an example application, we 
consider explicitly Higgs boson production in association with a single jet. This calculation 
is of direct phenomenological interest for understanding the properties of the new Higgs-like 
state observed at the LHC [TBI [19] . It extends previous work on understanding the effect of 
resummation on the Higgs plus zero-jet cross section [Qj [TOl [11]. We resum the logarithms 
\n{Q/p^^°) through the next-to-leading logarithmic (NLL) level, where Q ~ rriH ~ p'^ and 

is the transverse momentum of the observed jet. We demonstrate that the residual 
scale variation of the theoretical prediction is drastically reduced by the inclusion of the 
NLL resummation, and that the NLL-I-NLO result provides reliable predictions over a larger 
kinematic range. Since the factorization theorem we derive is valid for both more jets and 

^Wc note that different schemes for counting logarithms are employed in the literature; we specify in 
detail the order-counting scheme we use in Section [2j 
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other color-neutral particles, our result also serves as a framework for how to augment a host 
of fixed-order calculations with resummation of a class of large logarithms. 

In the context of recent work suggesting that it is difficult to extend resummation of 
jet- veto logarithms in the presence of the anti-Zcr algorithm beyond the NLL level [llj, a 
relevant question to consider is the required accuracy of the resummation to match experi- 
mental needs. Phrased more colloquially, how many N's of logarithmic accuracy are needed? 
In the example above, the numerical value of the leading logarithmic term in the experimen- 
tally interesting region is In {Q/Pt^°) ~ 2.5. This is not an overwhelmingly large correction. 
The ability to supplement fixed-order codes with simple analytical resummation of important 
sources of logarithms to the NLL level should be sufficient. This would allow these programs 
to be extended near regions of phase space where large logarithms appear, while still pro- 
ducing reliable central values and theoretical error estimates within a controlled theoretical 
approximation. We comment in more detail later in this manuscript on the numerical rele- 
vance of missing higher-order corrections that may arise at the NNLL level. We find that in 
the experimentally interesting parameter region for Higgs production, it is justified to resum 
the large logarithms associated with the jet veto. 

Our paper is organized as follows. In Section |2| we derive our factorization theorem 
using SCET. We apply our formalism to study Higgs plus one-jet production, and discuss 
our numerical results, in Section [3j We conclude and discuss future directions in Section |4| 
All formulae needed for resummation at the NLL level are given in the Appendix. 



2 Factorization and Resummation 



In this section we derive a factorization theorem for multi-jet production at the LHC in the 
presence of a jet veto. We discuss the resummation of the logarithms associated with the jet 
veto through NLL accuracy. We use pp — )■ Higgs + 1 jet via gluon-gluon fusion as an example 
to highlight the derivation procedure. The generalization to additional jets is straightforward, 
and is presented here. Our primary results are contained in Eqs. (14), (15), (26) and (27). 



2.1 Discussion of the jet constraints 

We focus on the case in which the jets are defined using the hadron collider anti-Zc-r algo- 
rithm [20]. The following distance metrics are used: 

Pi = VtI- (1) 

The anti-fcT algorithm merges particles i and j to form a new particle by adding their four 
momenta if pij is the smallest among all the metrics. Otherwise, i or j is promoted to 
a jet depending on whether pi or pj is smaller, and removed from the set of considered 
particles. This procedure is repeated until all particles are grouped into jets. We note that 
Ai?^- = Ar/?- + A0f -, where Ar^jj and A0jj are the rapidity and azimuthal angle difference 
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between particles i and j, respectively. R is the jet radius parameter and in practice is 
chosen to be around 0.4 — 0.5. 

After clustering the partons, we demand that the final state contains no additional jets 
with transverse momentum greater than a threshold p^^°. For Higgs searches, typically 
pjKto 25 — 30 GeV. Since p^^° is usually substantially lower than the partonic center of 
mass (A = p1^^"/\/§ <^ 1), the vetoed observables are usually very sensitive to soft and 
collinear emissions. 

Additional constraints beyond the jet veto can be imposed on the final state. In the 
following derivation, we require that the measured leading jets are all well-separated so that 
no additional small scales will be generated. We also assume that leading jet momentum 
p;f ~ rriH ~ and 1 ^ ^ \^ while ^log^i? <^ 1. The second of these requirements 
is necessary to insure that the measurement function factorizes into separate measurements 
in each of the collinear sectors. This is discussed in detail in the next section. The third 
requirement ensures that logarithms associated with the anti-Zcr parameter R need not be 
resummed. Given that p^^° ^ 25 — 30GeV and R ~ 0.4 — 0.5 for Higgs production,when the 
leading jet p;f ~ rnn ~ 125GeV, these assumptions are justified. 

2.2 Derivation of the factorization theorem 

We use SCET [131 ttH [13 [IE] to establish our factorization theorem. SCET makes 
manifest the infrared limits of QCD by re-formulating the QCD Lagrangian using soft and 
collinear modes whose momenta scale with a small power-counting parameter A in appropri- 
ate ways. For Higgs production, this parameter is of order pVf^° / for radiation outside the 
measured jet, and is of order R for radiation inside the measured jet. Consideration of the 
jet algorithm and jet veto lead to the following relevant degrees of freedom in the effective 
theory: 

• a collinear jet mode with momentum pj = '^nj + kj, where rij is the light-cone vector 
along the jet direction; 

• two collinear modes propagating along the beam axes a and b, with pi = + ki for 

i = a,b; 

• a soft mode with momentum kg. 

The residual momenta fcj, fcj and the soft momentum kg all scale as a/IA, while the large 
components of the three collinear momenta scale as uji ~ a/I. We note that no ultrasoft 
fields {k 

US 

are needed for the process we are considering here. Any phase space 
measurement Ai is assumed to be insensitive to these modes in the final state, so that 
^y_g M\Xus) {Xus\ = \Xus){Xus\ = 1 + C'(A), where X^s denotes the final-state particles 
with an ultrasoft momentum scaling. 

The leading-power SCET operator which mediates gluon-fusion Higgs plus one jet pro- 
duction is 

HO{x) = J2 e<"^'''^'^'''~'^''''^-''CXHSZ'sZs:i^^^^^ , (2) 
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where we have exphcitly written out the Lorentz indices and the color indices abc. H 
is the operator which creates a Higgs boson in the final state. C^g, is the hard Wilson 
coefficient which encodes the hard virtual ffuctuations. It can be obtained order-by-order 
by calculating the corresponding QCD diagrams. The n-collinear gauge invariant boson 
field HE] 

= ^ MiV^ + 9sA':^)Wn] 5^^p 5„,^ , (3) 
9s 

which creates (annihilates) a coUinear gluon in the final (initial) state, is built out of the 
coUinear Wilson line H] 



perm. 
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(4) 



V is the projective operator acting on the coUinear fields sitting to the right of it inside 
the parentheses. At leading power, the interactions between coUinear and soft fields can be 
eliminated through an operator redefinition B!^ — > S^"" B'^ , which results in the soft Wilson 
line Sf [15]. 

The cross section with a jet veto can be written as 

^ = / ^^e'''"-''{paPb\0\x)M\XaX,XjX:){X,XjX,Xa\OmPaPb). (5) 

spin X 

We have decomposed the final state into different sectors based on the momentum scaling. 
The operator O has been written in a factorizable form, but in order to proceed, we must 
demonstrate that the measurement operator M., which includes the jet clustering operation 
and jet vetoing, can also be factored into different sectors up to power suppressed correc- 
tions [21]. The factorizability of the measurement operator M. can be seen through power- 
counting the anti-ZcT algorithm metrics pij and pi. Recalling that the transverse momentum 
Pt for each sector scales as 

~ 0{l) , ~ ~ ~ 0{\) , (6) 

we have 

Pjj < Pj ~ 1 , pjs ~ , pja ~ pjb ~ log A"^ , 

Pss ~ Paa ~ Pbb ~ (A-R)""^ , Psa ~ Psb ~ Pab ~ (A-R)""^ log A"^ , 

Ps ^ Pa ^ Pb ^ A~^ . (7) 

These scalings indicate that for the jet-radius parameter R not too large {R ^ log A~^), the 
contributions from the mixing between the jet and beam sectors, and the mixing between 
the soft and beam sectors, are power suppressed. Jets tend to form separately within each 
sector. As long as i? ^ 1, radiation coUinear to the jet direction will be combined before it is 
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clustered with soft radiation. This means that the soft radiation is insensitive to the details 

of the coUinear radiation except for the jet direction. Therefore, A4 can also be factored 
between these two sectors. The measurement operator can therefore be factored as 

M = MjMsMaMb, (8) 

up to power-suppressed corrections in p^?*" and R. The individual Ma will be given in the 
Appendix. 

Plugging in the definition of the operator O, the cross section can be written as 



da 

d^H 8s 

spin ua 



I 

X 

Xa 



I 

I 

E(^^«l^nr (^n„, <, - (0) 

Xa 
I 

X E(p,|S£^^ (x-^, <, 0„J {X,\5{u, - P,)Sn„nX/^ (0) 

Xb 

I 

nj Xj 

Here, means summing over the final states with the restrictions imposed by M.. We 
have suppressed Ai for simplicity. In the last line, indicates the need to sum over 

all possible directions of different n.y-jet modes. We have removed several Kronecker-deltas 
using the discrete sums. The remaining ones have been turned into integrals J duA and delta 
functions (){uja — Va), after combining the residual momentum = ua ■ ^a with the label 
momentum uja in the coUinear sector A for A = a,h, J. We note that the coUinear sectors 
do not depend on x'^ = ua ■ x. 

We further modify this expression by performing a translational operation in each sector 
and inserting several residual momentum operators J dkS{k — k) to remove the explicit 
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momentum dependence on the final states. This gives us 



da 



d$ 



H 



^^J2Y1 J / / ^'^«/ ^^'^^^ 

spin HA 

I 

I 

Xa 
I 

Xb 
I 

c m'nTmxj){Xj\5{k] - kj)5'ikj - kj)5{coj - pj)B>;:f'^{om , (lo) 



where dk = dk^d'^k^. We now drop all residual momenta of order A or higher, keeping 
homogeneously only the leading power terms in the exponent (dropping the contributions 
of order A means that we ignore the recoil effect in the transverse plane). We perform the 
integration over the k~^ component in each coUinear sector to indicate that no restrictions 
are applied on this residual momentum |^ This leads to 



^ E y dwA J] ^"^kj (27r)^(5^ (^y + yn;, - y nj - qh^ 



da 
d^H 

spin nj 



xC^S^^'K' j d%{0\SZ'SZ'S:jiO)MsS\k,-k)S<'^'S^fS<f\0) 

X I d'ki{pa\B:':^{o)MaS\k^-k^)6{u,-v.)B:y^m^^^ 

X J d2A:,^(p.|SffK0)A^55'(A:,^-^,^)5(u;.-n)<'^'^(0)b) 

X (0|i3^f (0) MjS'ikj - kj)5{coj - Pj)B^J;<mo) . (11) 



To reach the formula above, we have summed over the final states using Y2x ~ 1 

and have re-inserted the measurement operators AiA- We next define the beam and the jet 
functions 

f±g/p{z, P^t"'", R)= j d^k^ Yl -^{p\BnO) Mb 6\k^ - k^)5{uj - P)S„,„,.(0) \p) , 

spin 

J{R)gT'6-' = 2i2nf{-ujmB^4{0)Mj6\kj - k j) 6 {uj -Vj)B(;' (0)10) , (12) 



^We note that this assumes that no pseudorapidity cut is imposed on the observed jet. It is straightforward 
to remove this constraint if desired; for simpHcity of presentation we do not do so here. 
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and the soft function 



s{nj,R,prn = J d%{o\sz's':,i^s:^;{o)Ms6\h - ^.)5r'^nf ^n;'(o)io) • (13) 

We emphasize that the beam and jet functions are well-defined only after the soft zero-bin 
subtraction has been properly performed [21]. 
Finally, we arrive at our result 

da = d$Hd$j J dxadxb— (^ p^2^_ {271)^5^ {qa + qb ~ qj - qu) 

XTt{H ■ S) f±g/pS^a,pT", R)f±g/p,{Xb,P^T"'", , (14) 

with the help of identifying J2n 2(2n)»uj J^"^^} massless particle phase space d$j = 



g|.2^dgjdn. We define the hard function H = CC"^ . The trace is over the color indices. 
Since p^^" is much larger than Aqcd, the beam function can be further matched onto parton 
distribution functions at the beam scale fiB ~ p^^°: fi^± = Xij ® fj{x). The matching 
coefficient Ijj can be calculated perturbatively. We have chosen our normalization so that at 
leading order, Xjj = 5(1 — x)6ij, J{R) = 1 and S = 5"" 6^^ 5'^^ . The factorization reproduces 
the tree-level gg Hg cross section. We note that for NLL resummation, only the leading- 
order matching coefficients are needed. In the Appendix, we will present the NLO results 
for the jet and beam functions, up to 0{R) corrections. 

The formalism is easily generalized to processes with an arbitrary number of jets and non- 
strongly interacting particles in the final state. All the arguments go through identically. 
We find 

da = d$H,d$j, J^i'^H,, '^jj Yl I da^adxb^ (27r)^5^ I + ^6 - gj, - J] ) 

a,b \ i c / 

_ _ n 

■ S)Ia,.j. ® /..(^a)X,,,,, ® f,,{Xb) n JjXR) , (15) 
spin color i 

where d^Hc and d$j^ denote the phase space measures for the color neutral particle and 
the massless jets Jj, respectively. J-'($//^, $j.) includes all additional phase-space cuts other 
than the pt veto acting on and the n hard jets, which should guarantee well separated 
n-jet final states (nj. ■ nj^ 3> A). The measured jet pj. should be much larger than p^^°. 

We note that possible issues arise when attempting to extend the resummation to the 
NNLL level, as pointed out in Ref. [11]. For R^ X, corrections of the form 

alR^lnX (16) 

appear, which prohibit soft-coUinear factorization. In the limit i? ~ A, clustering logarithms 
of the form 

a^lni^lnA (17) 
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which prohibit even NLL resummation occur. Let us study the numerical impact of these 
terms for the parameter values relevant for Higgs production. Setting R = 0.4, ttlh = 126 
GeV, and pf^° = 25 GeV, we find = 0.16, In (1/A) = 1.6, and In (l/R) = 0.92. For R = 
0.5, In (1/-R) = 0.69. It is clear that the corrections of Eq. (16) are indeed power-suppressed 
for experimentally-relevant value of i? ~ 0.4 — 0.5. They can be obtained to sufficient 
accuracy by matching to fixed-order results. There is also a hierarchy In (1/A) > ln(l/i?); 
the clustering logarithms are not large for relevant jet parameters. An eventual inclusion 
of the leading clustering effects by combining the resummation with a NNLO calculation of 
Higgs plus one-jet production should be sufficient. We therefore believe that it makes sense 
to study the resummation of jet-veto logarithms in Higgs production, and proceed with our 
analysis. 



2.3 NLL Resummation 



Each ingredient in Eqs. (14) and (15) describes fluctuations with a particular momentum 



scaling. When the hard, jet, beam and soft functions are calculated near their natural scales, 
no large logarithms will arise. The typical scales for each sector are 



Pt 



Hj ~ PtR , fJ^B ^ fJ'S 



veto 

Pt ■ 



However, calculating the cross section requires all the functions to be evaluated at the same 
factorization scale n, which generates large logarithms of the following ratios: 



Pt 
/i 



p^R 



^veto 

Pt 



(19) 



These large logarithms can be resummed by evolving each function to the scale /i using the 
renormalization group (RG) equation 



(20) 



The anomalous dimension 7^ can be extracted most easily from the e poles of each function 
calculated using dimensional regularization. 

However, like for small-g^ resummation, a subtlety arises because of the identical virtu- 
ality shared by the collinear and soft degrees of freedom [2H . Various efforts have been 
made in SCET to regulate this rapidity divergence [251 [26], |27] , which all have shown be able 
to reproduce correctly the NLO flxed-order QCD singularities. In our current approach, we 
adopt the formalism proposed in Ref. [25j. We regulate the extra divergence by modifying 
the collinear and the soft Wilson lines as follows: 



exp ( -gs 



perm. 



'V 



w 



perm. 



exp -gs 



'V 



w- 



n ■ Ar 



n ■ 



(21) 
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where the bookkeeping parameter w will be set to 1 as 77 —i- 0. The effective rapidity cut-off 
V and the new parameter i] play similar roles as and e in dimensional regularization. The 
corresponding rapidity-RG equation commutes with the conventional one for the beam and 
the soft functions: [/^ J;, '^J^] =0, 

= lu{iy)FBA^) ■ (22) 
The natural i/-scale choices for the beam and the soft functions are 

VB--Xa,b^S, Vs--pT°- (23) 

Since the physical cross section is /x and v independent, the anomalous dimensions must 
obey the consistency conditions 

7b + 75 = 0, (24) 

up to power corrections of order A and effects due to finite jet (beam) separations [22|. We 
will use these conditions to extract the anomalous dimension for the soft function. The 
general solution to the RG equation can be formally written as 

v) = U (/i, z/, /io, i^o)F{fio, t'o) • (25) 

The explicit form of the evolution kernel U for each function along with all details needed 
for NLL resummation are given in the Appendix. 

Recalling that for NLL resummation only the tree-level Wilson coefficients are needed, 
we find the following simplified expression for production of a color-neutral particle plus one 
jet: 



ab 



xUhAI^^ AIh) Us,kilJ', V, /i5, ^s) 1-B,a,b{f^i I^B, l^B, Xa, Xb) Ttjifl, ^j, R) . (26) 

For a more general process, the evolutions of the hard and the soft functions will usually 
induce operator mixing in color space. Therefore, the NLL cross section for multi-jet pro- 
duction reads 

dciNLL = ^ dXadXbTl HI''q^''''^''\iIh) S UH,{k}{fJ', fJ'H)Us,{k}{fJ','^, IJ'S,^s) 

fa{fJ'B,Xa)fbifJ'B,Xb)lB,a,bifJ',l^,fJ'B,l^B,Xa,Xb) ]^ 7^ (/i, /i J, -R) . (27) 



ab 
X 



We note that due to the separation between the scales inside and outside the allowed jets, 
a full refactorization or a refactorization ansatz of the soft function [221 123] may be helpful 
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in improving the resummation accuracy. Since this involves Ini? effects and is therefore 
moderate for the experimentally-interesting p^^" and R, its will be left for further studies. 
There also exist non-global logarithms [28] beginning at the NNLL level whose resummation 
is beyond the scope of our formalism presented in this work. 

Different schemes are used in the literature to determine the accuracy of resummation 
prescriptions. When calling our result NLL, we use the order-counting defined in Ref. [B]. 
Denoting the large logarithms associated with the veto scale generically as L, our NLL 
captures the two leading logarithms at each order in as- We correctly obtain a^L^ and asL, 
a^L^ and a^L^, a^L^ and a^L^, and so on. To obtain the next tower of logarithms (a^L^, 
a^L^, etc.), we would need to include the jet, beam and soft functions at 1-loop. This would 
correspond to NLL' in the language of Ref. [H]. All ingredients are currently known for this 
extension except for the one-loop soft function, which is not difficult to obtain. We plan 
to include this term in future detailed numerical studies. Control over the next tower of 
logarithms [a^L, a^L^, etc.) would require the two-loop non-cusp anomalous dimensions '-fx, 
and would correspond to NNLL in the notation of Ref. p]. 



3 Numerical results for Higgs+jet production 

We present in this section numerical results for pp — /i+jet at an 8TeV LHC with a jet veto 
imposed. We combine our resummation with the NLO cross section from MCFM to produce 
NLL+NLO results. In the numerics presented here we restrict the leading jet rapidity to 
\rij\ < 2.5 and veto all other jets with pT^i > p^*° and |77i| < oo for simplicity. Experiments 
typically only veto jets in the range \ri\ < 4.5. However, we expect that this boundary effect 
will be small |5|. It is simple to include this constraint if desired, as discussed in the previous 
section. We have included the gg, qg, and qq partonic scattering channels in obtaining these 
results. All relevant beam, jet and soft functions, as well as anomalous dimensions needed 
for this calculation, are presented in the Appendix. 

We begin by demonstrating that our formalism correctly sums all the next-to-leading 
logarithms of p^^" /pj- by comparing the expanded NLL production rate with the MCFM 
NLO result jT^. In the expanded NLL result, we have included the large non-logarithmic 
virtual corrections by using the full NLO hard function taken from Ref. [22] • The validity of 
our formalism is shown in Fig. [T| where we show the agreement between these two results in 
the region where \ogp^ / p^^" becomes large. This demonstration is based on the dominance 
of the log terms over the other contributions omitted in SCET in the small p^^" region. 



We next study the cross section for Higgs+jet production. Since Eq. (26) is only valid 
for pVf-^° <C p^, in order to give a prediction over the entire allowed kinematic range, we have 



to combine the resummed formula of Eq. (26) with the full NLO result. For this purpose 



we adopt the matching scheme proposed in Ref. [5] , in which the RG-improved cross section 
is taken as 

0- = I ) CXnlo(/W) - Z ((Texp(/U) - crLo(/u)) , (28) 
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Figure 1: Presented is the ratio of the expanded SCET cross section IcXexpl and the full NLO 
QCD calculation from MCFM ctmcfm, |c"exp — o"mcfm|/|o"mcfm|- We have required the leading 
jet pj' > 120GeV, |?7j| < 2.5, have set R = 0.1 and have made pl^^" as low as O.OlGeV. The 
excellent agreement between the expanded SCET prediction and the MCFM cross section 
for such low pl^^° implies that our formalism catches all the NLL structures. 



where Z = {\ — P^t^° I Pr^to) ■ "^lo is the LO cross section and ctnlo is the cross section 
calculated through NLO using MCFM. ciexp is obtained by expanding the resummed cross 
section ctnll in Eq. (26). We postpone rigorous study of the uncertainty induced by the 
choice of matching procedure to future work. In evaluating the resummed production rate 
c^NLL, we fix 



ix,v~sr-^^{x,v^sf-^''{p'^f-^' 



fij = PtR , 



(29) 



These choices minimize the logarithmic dependence in the hard, jet and the beam functions. 
When we use this set of scale choices, the cross section is also independent of the rapidity 
scale v. 

In Fig. [2| we demonstrate the improvement obtained with NLL resummation by showing 
the dependence of both the NLO and RG-improved integrated cross sections on different 
choices of p^^" . We have set p^ > 120GeV, fi = mn = 126GeV and have taken the anti-fcy 

120 GeV in Z. We use 
We can see that for small 



parameter R = 0.4. For the matching in Eq. (28), we use 
the MSTW2008 NLO PDF set ED] with two- 



,max 
veto 



oop as running. 

values of p^^°, the fixed-order cross section becomes negative while the NLO + NLL result 
remains positive. We also vary the scale /i from m///2 to 2mH to estimate the theoretical 
uncertainty. We can see that resumming the large logarithms greatly reduces the residual 
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Figure 2: Shown are the NLO+NLL prediction and the NLO cross section as a function of 
pJ^^° for p'^ > 120 GeV. The sohd-blue curve represents the RG-improved cross section in 
Eq. (28). The dashed-red hne shows the NLO result for the scale choice /i = mn, while the 
dotted red hues show the NLO result for = ttihI'^ and = 2mH ■ The blue band reflects 
the scale uncertainty of the RG-improved rate. The band boundaries are set by the values 
at the scale choices mn/'^ and 2m h- 



scale variation for the experimentally relevant values p^*" 25 — 30 GeV, leading to a more 
reliable prediction. The fixed-order cross section exhibits little scale variation for ?a 55 
GeV. Similar behavior is observed for the Higgs plus zero-jet cross section [31 H], and was 
argued to result from an accidental cancellation between several higher-order corrections 
with different origins. The same argument holds here for the Higgs plus one-jet result. As 
p1^^° becomes large, the fixed-order and resumed cross section coincide. Since the separation 
between the hard scale and p^^" vanishes in this limit, this behavior is expected. 

To gain some intuition regarding how low in pj. the resummation of jet-veto logarithms 
leads to a difference from fixed-order, we fix p1^^° = 25 GeV and integrate over the leading-jet 
transverse momentum subject to the constraint p;f > p^min- stress that some caution 
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must be exercised in using these results. Our effective theory framework is only vahd when 
the hierarchy p^^° ^ p;f exists. When p^^° ~ p^, our result reduces to the fixed-order 
result, which contains the ratio of scales *° <^ ttlh- A different effective theory framework 
consisting is needed to resume logarithms of this ratio. With these caveats stated, we plot 
in Fig. [3] the NLO and NLL+NLO results as a function of the minimum allowed jet pt 
for p^^° = 25 GeV. Significant differences between both the central values and residual 
scale uncertainties persist down to low values of p™", indicating the need to augment the 
fixed-order results with resummation over the entire pj^ region. 




Figure 3: The NLO+NLL prediction and the NLO cross section as a function of p^*" for 
pt^to _ 25GeV. The solid-blue curve represents the RG-improved cross section in Eq. (28). 
The dashed-red line shows the NLO result for the scale choice /i = ttlh, while the dotted red 
lines show the NLO result for /i = mH/2 and /i = 2m u . The blue band reflects the scale 
uncertainty of the RG-improved rate. The band boundaries are set by the values at the scale 
choices m///2 and 2m h- 
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4 Conclusions 



In this manuscript, we have estabhshed a formahsm for the production of color-singlet parti- 
cles produced in association with an exclusive number of jets at the LHC. Using effective field 
theory techniques, we have proven the factorization theorem of Eqs. (14) and (15), which 
allows us to resum large Sudakov logarithms of the form Xogp^^" /p^ to all orders. We have 
focused on Higgs production in association with a jet as an example. We have demonstrated 
by the excellent agreement between the expanded NLL result and NLO QCD cross section 
that our formalism correctly captures all relevant large logarithms in the p'^*'° — )■ limit. 
The scale uncertainty of the cross section is greatly reduced by inclusion of the resumma- 
tion. By matching our results with MCFM, we provide a NLL+NLO result for Higgs+jet 
production valid over the entire kinematic range. With our results, it is easy to supple- 
ment fixed-order results for vector boson or Higgs boson production with the resummation 
of jet-veto logarithms in order to provide predictions valid throughout phase space. 

Several future directions remain to be pursued. With our results it is possible to improve 
the theoretical predictions for Higgs plus one or two jets. This will be of great phenomeno- 
logical importance as the properties of the new state discovered at the LHC are further 
analyzed. We plan to further investigate these phenomenological application of our formal- 
ism. It is also interesting to study the effective theory valid when both the veto scale and 
the leading jet pt are smaller than irtH- Beyond NLL, extra clustering effects will enter the 
cross section, which are not completely understood [TT] yet. We have argued that the nu- 
merical impact of these terms should be subdominant to the effects studied here. However, 
a NNLO calculation of the jet, beam and soft functions would be helpful in determining 
whether this formalism can be extended beyond NLL. In the current work, we have simpli- 
fied the calculation of the anti-fcj- jet function by keeping the leading R contributions only 
in the resummation formulae. It will be phenomenologically interesting to include higher 
order corrections in R to improve the accuracy. We have also neglected for simplicity the 
boundary effects due to the experimentally finite i] range. These issues will be addressed in 
future detailed phenomenological studies. 
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A Fixed-order jet and beam functions 

In this Appendix, we list all ingredients needed for NLL resummation. We start with the 



NLO calculation of the jet and the beam functions, whose definitions can be found in Eq. ( 12 ). 
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The anti-kT jet function is calculated using the measurement function 



Mj = e(Ar;/, + A0f < R') + 0{prn • 



(30) 



We explicitly calculate that the jet collinear radiation leaking outside the jet is power- 
suppressed by p"^^" after correctly subtracting the soft zero-bin contributions. Since numer- 
ically R ^ 1, we can simplify the measure using 



Ar^J + = 2 cosh{Ar]ij) - 2 cos(A(/)y ) + 
In this limit, the NLO jet functions for gluons and quarks become 



(31) 



2n 
2n 



Ca 



67 _ 37r2 
13 Stt^ 



9 2 



+ C(i?' 



+ 31og-^ +2 log' 



+ 0{R^) . 



2 4 ' """"p^R ^ pJj,R 

The measure for the beam function with one emission is 

Mb = Q {kT, < pT°) (l^^l < ^cut) + > Vent) . 



(32) 



(33) 



Experimentally, rjcut ~ 4.5. For simplicity, we set r/cut = oo here. We note that this difference 
does not affect the anomalous dimension of the beam function, it only changes the finite 
part. The calculation is performed using the 't Hooft-Veltmann scheme. The NLO matching 
coefficient X for the beam function is found to be 



as{iJi)CA 
27r 

2^r 
as{iJ,)CF 
27r 

2tt 



4 log ^ log —5(1 -z)- 2p,g{z) log , 



4 log log - ^) - 2P..(^) log + (1 - ^) 

Pt 



-2pgq(^)l0g 



n-p 

^veto 

Pt 



Pt 



+ z 



-2P,.(^) log -^ + 2^(1-^) 
Pt 



(34) 



with 



2z 



{1-z). 
1 + z^ 

1 + (1-^)' 



+ 2^(1 - z) + 2 



1-z 



Pgg{z) ^l-2z + 2z^ . 



(35) 
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From these fixed order calculations, we can determine the anomalous dimensions used for 
RG evolution. The anomalous dimension for the jet function is given by 



7j. = 2r^^,pT- log 



a 



For the beam function, it is 



i'b 



2rcuspX!f, log , 



SFcuspT; lo^ 



n-p 



+ 1B, 



where = Ca for gluon and Tf = Cp for quark. Here, 



as (asY 



+ .... 



(36) 



(37) 



(38) 



B RG evolution 

The evolutions of the jet and the beam functions are given by 

-2TjArifij,fJ.) 

TZj. = exp[-2T^^S{fij,fi) - AjXf^j,fi) 



u 



B,a 



exp 



B. 



exp 



-T„Ar(/iB,/i) log — - AB^ (^B,/i) 



(39) 



We have defined lB,a,b = UB,aUB,b in Eqs. (26) and (27). For the Higgs production process 
considered in this work, we have, the following color identities: for the ggg channel, Tj ■ Tj = 
-Ca/2; for the gigac/s channel, Ti • T2 = -{Cp - Ca/2) and Ti ■ T3 = T2 ■ T3 = -Ca/2. The 
anomalous dimension for the hard function can be found in Ref. [21] • The solution to the 
RG equation is 



H 



exp 



2 T-S{fiH,f^) - 2AH{fiH, /i) + 2Ar(/iH, /i) Yl 



n 



where Ai?^- = 2 (cosh(A?7.y) — cos(A0jj)) for i.j 7^ a, 6, Ai?, 



Ai?: 



■ib 



(40) 



e^' and 



Ai?^^ = 1. Also = if i G J, otherwise Ua = Xa\fs- The soft anomalous dimension is 



determined by the consistency relation in Eq. (|24|), which gives the solution 

Us 



exp 



-2 Y T^Sifis, /i) - A,(^„ /i) - 2Ar(/i„ /i) ^ ^-^ log ARl 



ieB 



(41) 
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For the NLL resummation, we need 



(42) 



and 



ro r 47r 



4^2 {asii^i) 



1 - 



J - logr) + - I) (1 - r + logr) + A log^ ,| , 

(43) 



where r = as{iif)/ots{iii) and 



^0 = 



3 
34 



'1 — 



ri = 4 



^2 



3 
20 



20. 



3 

67 TT^ 
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(44) 

(45) 
(46) 

(47) 



Aj/B, Ah and are needed at leading order, and can be obtained by substituting the Fq 
in with the corresponding 70 and expanding in a^. For Aj/b we have 



1o' = QCf, 



and 



(48) 
(49) 



7o^ = 2^0 , 

for quark and gluon jets or beam functions, respectively. For = X^j7i, we have 

7o' = -3Cf, (50) 

and 

7o' = -A, (51) 

The two loop 75 can be extracted via the relation 2^h + 7Ji + 75 + 7b„ + Jb^ — 0. At one 
loop, 7^ = 0. 



C Expanded results 



Expanding out the resummed result will give us the fixed order singular term up to NLL. 
We find 



H _ 



9 * 9 



roE^iog^4 -roj:7;^iog^ + 7o^ + nA 
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where we have assumed at the leading order Hq oc a" and defined L = logi^jj,/ hh)"^ where 
is of order p^f , and Ai?f^- and uji have been defined previously. 
The expanded jet evolution gives 



Ji = l + 



47r 



roT^^iog^ 



J, 



To 
2 



(53) 



where L = log(/i//ij)^. We have checked that it reproduces singular pieces of the fixed-order 
NLO calculation of the jet function. To see this we note that the first two NLO terms can 
be combined to give log^ —j^ — - + (9(A°) as long as fij is of order p^'i?. 

The singular terms of the beam coefficient are given by 



1 + 



as(/i) 
An 



TT 



-Pa 



[z) log — , 



^oT^ log 



log ^ + log ^ log ^ 1 + 7o''" log 



FT 



fJ'B 



OJt. 



..Ba 



f^B 



(54) 



where Paa'{z) is the normal splitting function. This result again agrees with the fixed order 
calculation if we choose /is ~ p^^" and vb ~ i^a- 
Finally, the soft function is 



S 




47r 



^ 5: Tf log + ^ ri log ^ + r„ ^ 



aeB 



(55) 



where L = log(/i//i5')^ and fis ^ ^ p^^°. We note that once we combine the soft and the 
beam functions, the u dependence in these two contributions cancels, as required. If we set 
f^B = f^s = Pt^°, the Us and ub dependence also disappear. We note that numerically, the 
logi?^ X L term is of order 0.1 x L, which we interpret as a single- log term. 
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